Real-time simulation of finite frequency noise from a single electron emitter 
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We study the real-time emission of single electrons from a quantum dot coupled to a one dimen- 
sional conductor, using exact diagonalization on a discrete tight-binding chain. We show that from 
the calculation of the time-evolution of the one electron states, we have a simple access to all the 
relevant physical quantities in the system. In particular, we are able to compute accurately the 
finite frequency current autocorrelation noise. The method which we use is general and versatile, 
allowing to study the impact of many different parameters like the dot transparency or level posi- 
tion. Our results can be directly compared with existing experiments, and can also serve as a basis 
for future calculations including electronic interactions using the time dependent density-matrix 
renormalisation group and other techniques based on tight-binding models. 

PACS numbers: 73.23.-b, 73.63.-b, 72.70.+m, 
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I. INTRODUCTION 

The study of electronic transport in mesoscopic sys- 
tems - where small size, low temperature and careful 
fabrication of the sample ensure that quantum coherence 
is preserved - has shown tremendous progress in the past 
decades, providing a deep understanding of the behaviour 
of these systems^ Recent progress opened the way to 
the study of these systems on short time-scales, by look- 
ing at either the real-time dynamics or the high-frequency 
fluctuations^— A recent experiment has demonstrated 
the feasibility of single electron emission, where a single 
electron is periodically emitted from a quantum dot into 
an edge state of a 2d electron gas in the quantum Hall 
effect (see also Ref. [||). Measurements of the current 
and finite frequency current correlations in this system 
confirmed that the system indeed behaves as a single- 
electron emitter £ 

In this work, we study this system using real-time 
numerical simulations. Modeling the edge state and 
the quantum dot as a tight-binding chain without elec- 
tronic interactions (see Fig. [T]), we compute numerically 
the time evolution of the system, from that of the one- 
electron states, when a time-dependent gate voltage is 
applied to the dot. We are able to study all the aspects 
of single-electron emission, by calculating the average 
of all the relevant physical operators. In particular, we 
show that these real-time calculations allow us to com- 
pute accurately the finite frequency autocorrelation noise 
of the emitted current, which has been used as the def- 
inite experimental proof of single-electron emission. 9 In 
the optimal emission regime, our method provides very 
good agreement with the analytical results for the finite 
frequency noise obtained from a semi-classical models, 
which proves the power of real-time simulations in this 
context. More importantly, our approach can be used to 
study several different regimes where no analytical results 
are available. 

The goal of this work is twofold. First, as the single- 
electron emission experiments are currently performed 



in a regime where electronic interactions in the one- 
dimensional channel are small and can be neglected, the 
results obtained within our real-time simulation for non- 
interacting electrons are relevant as such, and can be di- 
rectly compared with experimental data. We stress how- 
ever that fermionic correlations are accounted for exactly 
with our approach. The method is quite versatile, and 
allows to easily study the effect of various parameters, 
like the transparencies, the temperature, the exact shape 
of the driving voltage, etc. Moreover, the real-time simu- 
lation gives access to appealing visualisation of the elec- 
tronic transport in these systems, with easy access to the 
real-time density in real and energy space, as well as the 
shape of emitted wavepackets. As a means to get results 
for the non-interacting system, our method thus appears 
as complementary to calculations based on the Floquet 
scattering theoryiir— , with equivalent results but differ- 
ent strengths and weaknesses. 

Second, this work is a necessary first step towards 
the use of more involved real-time simulations, such 
as time-dependent density matrix renormalization group 
(td-DMRG) and related techniques where electronic 
interactions can be taken into account, in such a time- 
dependent situation. 

The understanding of the role played by the different 
timescales, finite size effects and discretization, which our 
approach provides for non-interaction electrons, consti- 
tutes a prerequisite for the application of these more ad- 



FIG. 1: (color online) Sketch of the system: a tight-binding 
chain of length L, with a dot defined on one end (red sites, at 
the right of the td link) , on which a time-dependent potential 
V{t) is applied. 
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vanced and complex methods. In this regard, the present 
work is made in the same spirit as Ref. [24] (where shot- 
noise was computed from real-time dynamics) but actu- 
ally goes beyond in considering the finite-frequency noise 
in a time-dependent setup. 

This article is organized as follows. In Sec. [TT1 we de- 
tail the model which we have used, and we derive ex- 
pressions for the averages of physical quantities in terms 
of the numerical solution for the time evolution of the 
one-electron states. In sec. IIII1 we present our results, 
first considering the static properties of the dot (mean 
charge), and then extending these considerations to the 
time-dependent properties (emitted charge density along 
the chain, emitted current, and autocorrelation noise). 
Sec. HVl is devoted to discussions and conclusions. 



II. MODEL 

A. General formalism 

The system is modeled as a one dimensional (Id) tight- 
binding chain of length L, with the Hamiltonian: 

L 

H = - (cpi+i + h - c ) - *<* [Cjfid+l + h - c ) 

i>id 

where Cj and Cj are fermionic annihilation/creation op- 
erators at site j. The dot is composed of the sites j > jd, 
with a tunneling amplitude td to/from the dot, an on-site 
energy e^, and an applied time-dependent potential V(t). 
To perform periodic single-electron emission, V(t) is a 
periodic function of time (with period T), which in the 
optimal situation brings alternatively the highest occu- 
pied dot level above the Fermi energy in the leads (elec- 
tron emission), and the lowest unoccupied state below 
the Fermi energy (hole emission), see Fig. [21 In order to 
compute the evolution of the system, we consider it to 
be in equilibrium at t < with V(t) = const, while the 
periodic potential is applied for t > 0. 

For t < 0, the system is in equilibrium, and we can 
compute the L eigenstates of the Hamiltonian H, noted 
<j>i (i — 1 . . . L) with increasing energies E\ < E% < ■ ■ • < 
E^. The initial TV-particle state, at zero temperature, 
is simply obtained by filling the N lowest energy states. 
For non-zero temperature, the occupation of the states is 
given by the Fermi distribution. In the following, we will 
work with the chain at half- filling, with N = L/2. For 
t > 0, the evolution of each one-electron states 4>i(t) can 
be obtained by solving the time dependent Schrodinger 
equation: 

ih ^gl =H( t )m . ( 2 ) 



This can be achieved by numerical integration of the dif- 
ferential equation for an arbitrary potential V(t). 

As there are no electronic interaction in the system, 
each occupied electronic state evolves independently form 
the other states, and the TV-particle state at time t > is 
simply obtained by filling the same 1-electron state as for 
the initial state, but using now the time-dependent 4>i(t). 
Introducing 4>\{t) as the fermionic operator creating an 
electron in the single-electron state i at time t, the N- 
particle wavefunction of the chain at time t is: 

|*(t)> = ... 4 /2 (<) |0) 

= |1, 1,..., 1,0,0,. ..,0) t , (3) 

where |0) denotes the empty band. 

We want to stress that this method does not rely on 
any adiabatic approximation for the evolution due to 
V(t), and it can be applied with an arbitray potential 
V(t). In an nutshell, what we are doing is simply to take 
the system in a given initial state (t < 0), and to compute 
its time-evolution using the time-dependent Schrodinger 
equation. This can be done in principle for any system 
(e.g. even with electronic interactions), but would be im- 
practical at the level of the TV-electron state because of 
the huge size of the Hilbert space. But since we are con- 
sidering here non-interacting electrons, this can be done 
for each electron independenlty (Eq.©), thus relatively 
easily. 

For the sake of simplicity, we consider here the spe- 
cial case of a piecewise constant potential V(t), which 
allows to reduce the numerical integration of the differ- 
ential equation to simple matrix products (see below). 
Note however that the numerical cost for solving numer- 
ically the differential equation with a more complex V(t) 
is only slightly higher, so that in practice our method 
can be applied without difficulty to any reasonable form 
of V(t). We thus focus on a potential V(t) which at t > 
consists of perfect periodic steps (see fig.H]) and takes the 
form: 

!V- if t < 
V+ if + nT < t < | + nT (4) 
V- if ?-+nT <t<T + nT 

where n is a positive integer. In this case, knowing the 
eigenvalues and eigenstates of the system for V = V- and 
V = V+ is enough to get the time evolution of any one- 
particle state, by expressing it alternatively in the basis 
of these eigenstates. During a time interval where V(t) 
is constant, the time evolution is simply given by phase 
factors coming from the eigenenergies. Then the sudden 
switch of V(t) from V+ to V- (or vice-versa) is accounted 
for by a basis transformation from the V+ eigenstates 
to the V- eigenstates (or vice- versa), which ultimately 
amounts to a simple matrix product. By combining the 
trivial time-evolution for a constant V(t) with the change 
of basis when V(t) switches from one value to another, 
we easily get access to the one-particle state evolution for 
any time t. 
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FIG. 2: (color online) Top: Profile of the potential V(t) used 
in our calculations: V(t) is constant and equals V- for t < 
(system in equilibrium), then shows perfect periodic steps 
between V+ and V- for t > 0. Bottom: illustration of the 
electron emission process in the optimal regime, where V+ — 
V- = A (the level spacing of the dot), and the values of V+ 
and V- are chosen such that the Fermi energy cf of the chain 
is always in the middle of two energy levels. The left drawing 
shows the electron emission, when the driving potential brings 
an occupied level (black) above the Fermi energy; the right 
drawing shows the hole emission, when the empty level (black) 
is put below the Fermi energy by the driving potential. 



B. Calculation of operators average 

The fcrmionic operator for the one-particle state 4>i(t) 
can be written in full generality in the basis of the posi- 
tion states along the chain: 



i(4) = J2a hJ (t)cj, 



(5) 



where ctij (t) is a time-dependent unitary matrix. This 
matrix is the natural object coming out of the numerical 
evolution of the system, as it expresses the one-particle 
states at time 4 in terms of the position states. In par- 
ticular, the lines of the matrix aij(t) at 4 = corre- 
spond to the eigenstates of the Hamiltonian (JXJ) , which 
are naturally obtained in the position state basis, as the 
Hamiltonian itself. 

For convenience, we compute all the operator averages 
in the Heisenberg picture. The time-dependent position 
operators along the chain can then be written as: 



1. Average current 

The current operator on link k (between sites k and 
k + 1) at time 4 is: 



(7) 



4(4) = -i (4(*)cfc+i(i)-h.c.) 
and its average is: 

ih(t)) = -*(^Q!m,fc(*X,fc+l(*) 
n 

'¥(0) ^(0)^,(0) *(0)\-h.c.;) (8) 



with |*(0)) = |*> the half filled Fermi sea at time t = 0. 
Using Eq. ([3]), at zero temperature, the average of the <j> 
operators product is simply <5 mj „ for the occupied states 
(to < L/2), and for the empty ones (m > L/2), leading 
to: 

{/*(*)) =2 M«»,*(*Kfc+i(*))- ( 9 ) 

n<L/2 

At non-zero temperature, one has {(jy n (j) n ) = f(E n ), 
where f(E) is the Fermi function at energy E, and the 
average current becomes: 

L 

(I k (t))=2j2f(E n )Im(a n , k (t)a* ntk+1 (t)) . (10) 



2. Current correlations 

The correlation between the current at link k and time 
t\ with the one at link I and time ti is given by: 

S k Ah,t 2 ) = <*|/ fc (ti)/j(t2)|*> 

-<*|/ fc (ti)|*)(*|/i(*a)|*). (11) 

Proceeding along the same lines as for the current av- 
erage, and using Wick's theorem for averages involving 
four <j) operators, we get at finite temperature: 

L L 

s k ,i(h,t2) = (-1) E - 

m— 1 n— 1 

x (a m ,fe(ti)Q^ )fc+ i(ti) - a m ,fe+i(ti)a* )fc (<i)) 

x (a n .i(t 2 )a* m l+1 (t 2 ) - a rh i +1 (t 2 )a* m l {t 2 )) . (12) 



c k (t) = ^ 0^^(4)^(0) 



(6) 



where the one-electron states operators are taken at 4 
0. 



3. Mean charge on the dot 
The operator Q(t) describing the charge on the dot is: 
Q(t) = £ c](t)c 3 (t), (13) 

3>jd 



4 



leading to the average charge: 



j>jd m,n 

x(*(0)|^(0)^ n (0)||*(0); 

L 

i>id n=l 

The results from Eqs. (|T0|). (fT2|) and CEU) show that 
all the averages can easily be computed from the time- 
dependent matrix cti,j(t), which is directly obtained from 
the numerical computation of the time evolution. 



III. RESULTS 
A. Parameters and operating regime 

We present below the results which we have obtained 
to characterize the emission of single electrons in the one- 
dimensional chain. One must be aware that our approach 
uses a non-chiral system, where excitations can propagate 
to the left and to the right, in order to simulate a chiral 
system (an edge state of the quantum Hall effect). In our 
setup, the left-going excitations correspond to the ones 
coming out of the dot in the real system, while the right- 
going excitations are incoming onto the dot, and thus 
need to be avoided. Multiple reflections inside the simu- 
lated dot are equivalent in the real system to performing 
multiple round trips inside the dot. 

In the remainder of the text, we work with a dot com- 
posed of 5 sites, located at the end of a much longer chain 
(~ a few hundred sites, typically 500). This size of the 
dot is the result of a compromise. On the one hand, the 
dot has to be as small as possible compared to the chain, 
in order to minimize finite-size effects. Indeed the level 
spacing in the chain must be negligible compared to the 
level spacing in the dot, to appear as a continuum. On 
the other hand, since we apply a time-dependent voltage 
of the order of the dot level spacing, the latter needs to 
be small compared to the bandwidth of the tight-binding 
chain, in order to minimize the effects of the non-linear 
dispersion relation of the chain. As the level spacing is in- 
versely proportional to the size of the dot, one thus needs 
to consider a large enough dot to meet this requirement. 
In practice, we found that using a dot of 5 sites was a 
good compromise, in particular because we were able to 
perform the current measurement near the dot, thus re- 
ducing the spreading of wavepackets during propagation, 
which arises from the non-linear dispersion. 

We consider the chain to be at half-filling, and work 
with physical dimensions corresponding to h = e = 1. 
The tunneling amplitude along the chain (first term in 
the Hamiltonian ([T])) has been taken as the unit of energy, 
and as a consequence the Fermi velocity at half-filling is 
equal to 2 sites per unit of time. This means that an 



excitation takes a time L to reflect at the boundary and 
come back to its starting point along a chain of L sites. 
In particular, the current measured at a given link due to 
the passage of an excitation will get completely spoiled 
after a time L by the reflection of that same excitation 
at the end of the chain. One can thus expect that the 
numerical quantities which we compute will be reliable 
up to a time ~ L only. 

Unless specifically stated, all the results presented be- 
low were obtained in the optimal regime for electron 
emission. In this regime, the amplitude of the driving 
potential, V+ — V—, is taken to be equal to the dot level 
spacing A. When the driving potential switches to the 
value V+, one energy level of the dot is put at an energy 
ep + A/2 (where ep is the Fermi energy of the chain), and 
an electron is emitted at this very energy (see the bottom 
left panel of Fig. [5]). When V(t) switches back to V_, this 
same level is brought down to the energy ep — A/2, and 
a hole is emitted, i.e. an electron tunnels back to the dot 
level (see the bottom right panel of Fig. [5J . This parame- 
ter regime is the best suited to emit a true single electron 
and a single hole during each period, as has been shown 
experimentally 7 , and it is thus the most relevant one for 
our study. Note however that there is no difficulty in 
exploring other operating regimes with our method (for 
example, a dot level in resonance with ep). 



B. Static properties of the dot 

Before studying the properties of the electron emission 
by the dot, it is necessary to characterize the properties 
of the dot itself without a time-dependent potential V(i) 
applied. Fig. [3] shows the mean charge on the dot as a 
function of the on-site energy e^, for several values of the 
tunneling amplitudes to the dot t& at zero temperature 
(top plot), and for increasing temperature with t& = 0.2 
(bottom plot). For a small value of td at zero tempera- 
ture, we observe that the mean charge is quantized, with 
well-defined plateaus at integer value. This simply re- 
flects that when the dot is weakly coupled to the chain, 
the number of electrons on the dot is an integer: it can 
vary between 5 and 0, leading to five energy levels. In- 
deed, the density of states (not shown) which corresponds 
to d(Q)/dep = —d(Q)/ded has five narrow peaks at the 
position of the steps in (Q), corresponding to the five 
energy levels of the dot. As td is increased, the steps 
between the plateaus are smoothed due to the increased 
fluctuations between the dot and the chain, leading to 
broadened energy levels. Increasing the temperature for 
fixed td (bottom plot) has a similar effect: it broadens 
the energy levels, leading to smoothed steps. 

For the parameters which we have chosen, the level 
spacing of the dot A is of the order of 1. Yet, because 
of finite size effects, the level spacing A is not constant 
and depends on the levels considered. We thus focus on 
the electron/hole emission from the central level of the 
dot (ed — 0), for which the optimal emission regime is 
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FIG. 3: (color online) Mean charge on the dot as a function of 
the on-site energy of the dot e^, for several tunneling ampli- 
tudes to the dot td at zero temperature (top), and for several 
temperatures at td = 0.2 (bottom). For small td and tem- 
perature, the mean charge tends to be quantized at integer 
values, with sharp steps between the plateaus. When td or 
the temperature are increased, the steps are smoothed due to 
the increased fluctuations between the dot and the lead. 



obtained by taking V + = —VI 



0.6. 



C. Time-dependent charge on the dot 

We consider now the time-dependent charge on the 
dot when V{t) is applied, in the optimal emission regime 
{td = 0, V + = — V_ = 0.6). Fig0] shows the mean charge 
for this choice of V+, V- as a function of time (0 < t < T), 
for several values of the tunneling amplitude to the dot td, 
at zero temperature. Several observations can be made 
from these plots. 

First, when the potential V is switched from V- to V + 
(at t = 0), one electron leaves the dot, and the charge, up 
to a very good approximation, decreases exponentially: 



(Q)(t) ~ 2 + exp(-t/T) 



(15) 



with a characteristic time r which increases as td de- 
creases. Calculations in the continuous limit for a perfect 
dot predict: 7 



2tt 1 



(16) 



where D is the transparency of the dot. For the tight- 
binding chain which we consider, the relation between 
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FIG. 4: (color online) Mean charge (Q) on the dot as a func- 
tion of time when V(t) is applied for different tunneling to 
the dot td at zero temperature. This data is obtained for a 
chain of length L = 500, with a period T = 200. 



the tunneling amplitude td and the transparency is, 
the linear regime of voltage^ - — 



D(t d ) 



«3 



(17) 



Using this expression of D(td) and performing a fit of 
the data gives the function r ~ 6.1/D(td), with a rea- 
sonable agreement with the value of the level spacing in 
our system (6.1 ~ 2ir/A for A ~ 1). 

When the potential V(t) switches back to V- (for t > 
T/2), the charge starts to increase (as an electron comes 
back inside the dot, or equivalently a hole comes out 
of the dot), and the behaviour is symmetrically reversed, 
with a similar exponential growth back towards the initial 
charge. 

Second, the qualitative behaviour of the system 
changes as the ratio of the decay time r to the half- 
period T/2 changes. When r <C T/2, the electron is 
emitted/absorbed from the dot much faster than the half- 
period, and the emitted/absorbed charge is close to 1. 
This is the case on the two upper plots of FigJH where 
it is clear that the charge (Q) goes all the way from 3 to 
2 on the first half-period, then back from 2 to 3 on the 
second half-period. The electron is thus emitted or ab- 
sorbed with a probability close to 1 at each half-period. 
On the other hand, when r > T/2, the electron emis- 
sion/absorption is not complete during one half-period, 
which means that the probability that an electron is emit- 
ted or absorbed on each half-period is smaller than one. 
This is shown on the bottom plots of Figj4] 

Note that, in this regime (r > T/2), the variation of 
the charge during one half-period is smaller than 1. One 
has then to be careful with the choice of initial conditions. 
Indeed, we are trying to model a periodic system using 
a finite one starting from an equilibrium/ground state 
configuration. Depending on the initial parameters, the 
system will take one or several periods to relax towards 
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FIG. 5: (color online) Mean charge (Q) on the dot as a 
function of time when V(t) is applied, in the regime where 
t 2> T/2 (small tunneling amplitude to the dot). The left 
plot, obtained with a period T = 50 on a chain of length 
L — 1000 shows that starting at t = with (Q) = 3 does not 
correspond to the long time limit: after several periods one 
reaches a long-time regime where (Q) oscillates between the 
values 2 + Q+ and 2 + Q_. The right plot, obtained with a 
period T = 200 on a chain of length L — 500, shows that it is 
possible to access the long-time limit from t — by switching 
on V(t) at an earlier time ti„ < to get the correct charge 
{Q)(t = 0) = 2 + Q+ (see text for details). 



the permanent periodic regime (which is independent of 
the initial conditions). This is illustrated on the left plot 
of Fig[5j which shows (Q)(t), with r 3> T/2, for several 
periods (we chose T = 50 for a chain of length L = 1000, 
to be able to compute the evolution on many periods 
without suffering from reflections at the boundary of the 
chain) . It is clear that the extrema of the oscillations of 
(Q) tend to relax towards long term values 2+Q + and 2+ 
Q-. As it is numerically very inefficient to compute the 
evolution over many periods (it requires very long chains 
to avoid boundary effects), we have to modify the initial 
conditions in order to reach the periodic regime as soon 
as the time-dependent potential is applied. Knowing the 
value of the escape time r from the dot for our choice of 
parameters, and introducing s = exp(— ^-), the values of 
Q + and Q_ can be predicted analytically, by solving the 
equations expressing that the charge at the end of one 
period (at t = nT + T — 0) has to be the same that the 
charge at the beginning of the period (at t = nT + 0): 



D. Visualization of the time-dependent density 

The numerical results allow to obtain a visual represen- 
tation of the emitted wavepackets propagating along the 
chain. The left part of Fig. |5]shows an example of results, 
obtained for a chain of length L = 500, and an interme- 
diate value for the dot tunneling amplitude, tj = 0.4, at 
zero temperature. The chain is initially in equilibrium at 
half-filling (occupation =0.5 on all sites, except on the 
dot at the extreme right of the chain which contains 3 
electrons). Starting at t = 0, the electron is emitted, and 
appears as a left-propagating wavepacket, which initially 
has an exponential profile. As it propagates along the 
chain, this wavepacket broadens due to the non-linear 
spectrum of the chain. Ati = 100 = T/2 the emission of 
a hole starts, with similar properties as the electron. 

Performing a Fourier transform for each of the one- 
electron states gives us access to the density in the wave- 
vector space (or k space) as a function of time, which is 
shown on the right part of Fig. [51 Initially, the system 
in equilibrium, at zero temperature, has all the states 
for \k\ < k max occupied (occupation 1, the "Fermi sea"), 
while all the other states are empty. For t > 0, the elec- 
tron emission shows up as a well-defined peak above the 
Fermi sea. Similarly, the emission of the hole for t > 100 
appears as a dip in the Fermi sea. Both structures man- 
ifest at the same distance from hp, confirming that the 
electron and hole are emitted with energies symmetric 
with respect to the Fermi level. 

Looking at the density in fc-space for different tunnel- 
ing amplitudes td gives us information about the proper- 
ties of the emitted electron, and about the "quality" of 
the electron emission. Fig. [7] shows the electron peak for 
a small tunneling amplitude td — 0.2 and a larger one 
td = 0.6 (and other parameters as in Fig. [5]). A small 
transparency leads to a narrow peak which is well sep- 
arated from the Fermi sea, while a larger transparency 
creates a broader peak which tends to merge with the 
Fermi sea»ii 



E. Time-dependent current 



Q_ = sQ+ , Q+ = s Q_ + (1 - s) (18) 



giving Q + = 1/(1 + s), Q- = s/(l + s). In practice, in 
order to have the correct value of the charge (2 + Q+) at 
t = 0, we switch the potential V(t) from V- to V+ at a 
time ti n < 0, such that the value of (Q) at time t = is 
precisely 2+Q + . We are thus letting a fraction of electron 
escape from the dot for t < 0, in order for the charge to 
take the right value at t = 0. The result of this procedure 
is shown on the right plot of Fig{5j One can see that the 
charge has indeed the predicted periodic behaviour, with 
increase/decrease between the values 2 + Q + and 2 + Q- 
shown as dashed lines. 



The current coming out of the dot contains similar 
information to the time-dependent charge on the dot, 
as it is given by the time derivative of this charge. We 
recall that the current is computed on a link very close 
to the dot (a few sites away), in order to minimize the 
effect of the non-linear spectrum of the chain, which tends 
to spread the wavepacket propagating along the chain 
during its travel. 

FigJS] shows the current for the same parameters as 
Fig HI computed from Eq. ^ between sites 490 and 
491 for a chain of length 500. For t d = 0.5 and 0.3, 
the current shows a dip followed by a peak, due to the 
passage of the electron then the hole. The peaks have 
an exponential profile with some oscillations which cor- 
respond to the ones present in the charge (Q)(t) (we 
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FIG. 6: (color online) (Left) Density along the chain (mean occupation of each site) at different times for a chain of length 
L = 500, with a dot tunneling amplitude td = 0.4, at zero temperature, over one period of length T = 200. One clearly sees 
the emission of the electron at t = with an exponential profile, and of a hole at t = T/2 = 100. Because of the non-linear 
spectrum of the chain, wavepackets broaden during propagation. (Right) Same evolution shown for the occupation in k space, 
for negative k (left-going excitation). The emission of the electron shows up as a well defined peak above the occupied states. 




FIG. 7: (color online) Effect of the dot tunneling amplitude 
td on the electron peak in momentum space, for td = 0.2 
(left) and td = 0.6 (right), for a chain of 500 sites, with other 
parameters as in Fig. [6] A smaller transparency leads to a 
narrower peak well-separated from the occupied states, while 
a larger transparency gives a broader peak which tends to 
merge with the occupied states. 
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could check explicitly that computing d(Q)(t)/dt gives 
the same peaks). As td is reduced, the width of the elec- 
tron and hole wavepackets increases, and becomes larger 
than the half period. This is apparent on the plots for 
td = 0.2 and td — 0.1, where the current does not reach 
zero at t = 100 = T/2. There, when the external drive 
switches from V+ to V-, and the hole emission begins, 
the sign of the current changes abruptly. Bear in mind 
that we had to start the emission at time U n < (see 
Sec. IIII C[) in order to be immediately in the periodic 
regime at t = 0. This has the undesired consequence 
to reduce the oscillations of the current in the first half- 
period compared to the following ones, as can be seen 
on the bottom right plot of Fig. [5] (the current for later 
periods -not shown- is similar to the one of the second 
half-period). This is because, when using a non zero 



FIG. 8: (color online) Mean current (I) between sites 490 and 
491 as a function of time when V(t) is applied for different 
tunneling to the dot td at zero temperature. Data obtained 
for a chain of length L — 500, with a period T = 200. 



the charge has the correct "periodic" value at t = 0, but 
there is no abrupt change of the potential at that time 
unlike what happens at t = T/2,T, .... The values for 
the current during the first half-period are thus some- 
what inaccurate. For this reason, in the calculation of 
the current autocorrelations (see below), we use the cur- 
rent starting from the second half-period only. 
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F. Current autocorrelations 

The current correlations at finite frequency contain 
more information about the single emission process than 
the mean current or the mean charge do. We have com- 
puted the autocorrelation for currents on the same link: 

S k , k (ti,t 2 ) = (*|7 fc (ti)/fc(t 2 )|*) 

- <*|Ifc(ti)|*> (*|I fc (t 2 )|*) (19) 

for a link k close to the dot. Because the applied voltage 
is periodic in time, Sk,k(ti,t2) is a periodic function of 
the mean time (t\ + *2)/2. By taking the average over 
one period T, we obtain the autocorrelation noise which 
depends on the time difference only: 

S Kk (t) = S k , k {t' 1 t> + tf . (20) 

Performing a Fourier transform of Sk &(£) gives access 
to the autocorrelation noise as a function of frequency. 
This quantity has been measured experimentally for a 
frequency of the order of the frequency of the applied 
voltage V(t), as reported in Ref.|9j. Using a semi-classical 
model, valid in the optimal emission regime where the 
Fermi energy of the chain is halfway between two dot 
levels, Albert et al.— have provided an analytical for- 
mula for the autocorrelation noise as a function of the 
frequency to and the escape time r of the electron from 
the dot: 

^H = |tanh(|) T ^, (21) 

where T is the period of the applied voltage. A good 
agreement between this formula and the experimental 
results was obtained.— 

Fig. IH1 shows the results obtained for Sk,k(w = 27r/T) 
in units of the driving frequency (1/T), as a function 
of the escape time of the dot scaled by the half period 
t/(T/2), for a chain of length L = 600 with a period 
T = 200. The current correlations were computed at 
zero temperature, close to the dot, on the link between 
sites 590 and 591 (while the dot occupies sites 596 to 
600). Note that in this regime, the results are nearly 
insensitive to temperature, as long as it is much smaller 
than the level spacing in the dot. The first observation 
to be made is that the numerical results (red points) are 
very close to the predicted analytical formula (dashed 
black curve): the numerical curve is just slightly above 
the theoretical one (by s» 4%). 

The curve has a maximum for r ~ T/2, delimitating 
two regions which correspond to two different regimes. 
On the left of the maximum, when r <C T/2, the emission 
time of the electron or hole is much shorter than the half- 
period, so it is emitted with probability 1, and the noise 
comes from the uncertainty in the emission time, which 
decreases with r. This is the jitter regime, or phase noise 
regime. On the right, when r 3> T/2, the emission time 



of the electron or hole is larger that the half-period, so 
the emission probability is smaller than 1, which is the 
cause of the noise. This is thus the shot noise regime, 
or charge noise regime. Note that the deviations of our 
numerical results with respect to the theoretical curve are 
similar in both regimes. This is a good indication of the 
reliability of our calculations, in particular that the use 
of an initial time f j„ < 0, which is only needed in the shot 
noise regime (see Sec. IIIIC|) gives correct results. Larger 
deviations appear only for very small values of r/(T/2), 
which correspond to the large transparency regime when 
the escape time from the dot gets close to the time of 
a round trip inside the dot. The semi-classical model 
cannot be used in that regime and the analytical formula 
(ET1) is no longer valid. 

We have chosen the period T = 200 as large as possi- 
ble for the chain length L — 600, taking into account the 
finite size constraints. Indeed, the maximum time before 
an emitted wavepacket is reflected at the boundary of 
the chain and comes back spoiling the measurement of 
the current, is approximatively t max ~ 600 (because the 
Fermi velocity is 2 sites per unit of time). The time av- 
eraging procedure (Eq. ([2U1) ) requires at least one period 
T, while the integration for the Fourier transform needs 
at least another period T. Moreover, we had to discard 
the current obtained during the first half-period, as some 
of its features are different from those of the real periodic 
current (see the lower right panel of Fig. [Sj and discus- 
sion at the end of Sec. KITE]). Taking T = t max /3 ~ T/3 
is thus close to the maximum value which can reasonably 
be taken for the period T. Another constraint prevents us 
from taking arbitrarily large values for the escape time. 
The magnitude of the initial emission time U n < that 
we need to use increases with r, making it difficult to 
reach numerically large values of r, as it imposes to have 
a larger maximum time t max , and thus a longer chain. 
This is the reason why the largest value of t/ (T/2) on 
Fig. O is approximatively 2. 

In order to investigate the effects of the finite size on 
the difference between the numerical results and the the- 
oretical curve, we have computed the current autocorre- 
lation for increasing sizes, L — 300, 400 and 500, taking 
T = L/3 in each case. As shown on Fig. [10] (left), tak- 
ing a smaller chain, and thus a smaller period T, gives a 
higher current autocorrelation. We also see that the devi- 
ations from the analytical model for large transparencies 
(bottom left part of the curves) appears for larger val- 
ues of t/(T/2) for decreasing chain lengths L. This is 
because a given value of r/ (T/2) corresponds to a larger 
transparency for a smaller chain, which has a smaller T. 

The right plot of Fig. [10] shows the relative difference 
between the computed autocorrelation and the theoret- 
ical prediction, AS/S = {S(k^{2ir/T) - S fheor )/ S t heor, 
as a function of the escape time r. We see that this 
relative difference decreases from the 8-10 % range for 
a chain with L = 300 and T = 100 to approximatively 
4% for the L = 600 chain with T = 200. One can thus 
expect that using a longer chain with a larger T should 
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FIG. 9: (color online) Autocorrelation noise S of the current 
on link 590, taken at the frequency 2ir/T of the applied poten- 
tial V(t), as a function of the escape time from the dot. The 
noise is expressed in unit of the drive frequency (1/T), while 
the escape time is in unit of the half-period (T/2). The black 
dashed curve shows the theoretical prediction of the semi- 
classical model of Ref. [Tol . The red points show the results 
obtained with our calculations, for a chain of length L — 600 
and a period T = 200. 




FIG. 10: (color online) (Left) Same as Fig.[9]for chain lengths 
L — 300, 400, 500 and 600 (blue, orange, green and red, or 
from top to bottom). In each case the period has been cho- 
sen as T = L/3. Smaller chain gives a higher autocorrelation 
noise. (Right) Relative difference AS/S between the com- 
puted autocorrelation and the theoretical formula, for the 
same chain lengths, as a function of the escape time from 
the dot. The relative difference decreases when the chain 
length (and thus the period T) increases; for the longer chain 
(L = 600), the relative difference is about 4% 



bring the numerical results closer to the theoretical pre- 
diction, although a systematic difference, independent of 
the chain length, cannot be excluded from our results. 

Our simulations allow us to compute easily the auto- 
correlation noise for different operating regimes, which 
cannot be described by the semi-classical model. As an 
example, we have considered the resonant regime, where 
for the two values of V(t) there is a dot level resonant 
with the Fermi energy (see left panel of Fig. [TTj) . This 
regime is not expected to provide good quality emission 
of single electron, as two levels are contributing equally, 
at the same time, to the electron emission (the level at 
the Fermi energy, and the level above it). This is clearly 
visible in the computed autocorrelation noise (right panel 
of Fig. [TTJ) , which is much larger than the noise in the 
optimal regime. This regime is also much more sensi- 
tive to the temperature since a dot level is always at the 
Fermi energy: for small emission time r, the noise ob- 
tained at a temperature of 0.1 is larger than the one at 
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FIG. 11: (color online) Autocorrelation noise in a same setup 
as in Fig. in the resonant regime: dot levels are put at the 
Fermi energy by the external potential (see left panel). The 
curves with dots show the results obtained at temperature 
(red) and 0.1 (blue). The black dashed curve is the analytical 
prediction in the optimal regime. The noise in the resonant 
regime is much higher than in the optimal regime, showing 
that the resonant regime is not suitable for single electron 
emission. 



zero-temperature. 



IV. CONCLUSIONS 

In this work, we have shown that real-time simulations 
on a tight-binding chain are a powerful and versatile tool 
to study the single-electron emitter, giving the ability 
to compute accurately elaborate quantities like the finite 
frequency current-autocorrelation noise. Using a non- 
interacting system, we have shown that the time evolu- 
tion can be obtained simply and efficiently from the evo- 
lution of the one electron states, and we have derived the 
formulas for the useful physical quantities (charge on the 
dot, emitted current, current autocorrelations) in terms 
of the numerical object describing the one electron states 
evolution. This allowed us to study the characteristics of 
the single-electron emitter as a function of different pa- 
rameters (dot transparency, period of the AC drive, tem- 
perature, ...), giving us access at the same time to illus- 
trative visualizations (current, time-dependent density 
along the chain, etc.) and providing information about 
complex quantities. Our calculations of finite-frequency 
current-autocorrelation noise are in good agreement with 
theoretical predictions in the regime when these are ap- 
plicable (and thus also with experimental observations), 
and can be used to explore other regimes. We have de- 
tailed the limitations due to the finite size of the system. 

This work can be extended in several directions. First, 
for non-interacting electrons, we have only explored a 
small part of the parameter space and of the possible 
manipulations of the system. We have mainly focused 
in this work on the optimal regime for single-electron 
emission as well as the less effective resonant regime, but 
the method can also be used to study in detail other 
regimes. We have also limited our study to the case of 
a perfect piecewise constant time-dependent potential. 
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Our approach can easily be used to study the effect of the 
experimental limitations on this shape of potential, or to 
study different types of potential (for example V(t) ~ 
sin(i)) which would lead to the emission of electron with 
different energy distribution. It also has the ability to 
simulate a true "1 shot" experiment, where a single pulse 
(rather than a periodic drive) is applied on the dot to 
emit one electron, starting from an equilibrium situation. 

Our method can also be extended to study more elab- 
orate setups based on the manipulation of the emitted 
single electron. For example, a weak link in the middle 
of the chain, which would only partially transmit an in- 
coming wavepacket, would be equivalent to a quantum 
point contact. This could be used to perform a real-time 
study of tomography protocol proposals,— or to consider 
more complex setups like interferometers^ Placing dots 
at both end of the chain could also be used to study 
two-electron interferences processesj 29 ' 30 

Secondly, this work can be seen as a first step towards 
the study of the single electron emitter taking into ac- 
count electronic interactions. As it shows that using a 
discrete tight-binding model allows to get a faithful and 
accurate physical description of the single electron emit- 
ter, we can consider extending our calculations to more 
powerful methods like the time-dependent density-matrix 
renormalization group (td-DMRG) 16-23 which can in- 



clude electronic interactions. This work gives the nec- 
essary understanding of the role and limitations of the 
different parameters in the system (chain size, dot size, 
period of the AC drive, etc). Using td-DMRG (or similar 
techniques), we could study the impact on the emitted 
electron of strong Coulomb interactions inside the dot. 
We could also consider electronic density-density inter- 
actions between two edge states, which should be experi- 
mentally relevant as the experiments have been achieved 
so far with edge states of the Quantum Hall effect with 
a filling factor v = 2, where two edge states are propa- 
gating. 

Note added. A preprint (Ref. 15?) has recently ap- 
peared, where details of the Floquet scattering theory 
for the same setup are presented, together with the semi- 
classical model. 
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